Spin polaron in the J\ — J 2 Heisenberg model 
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We have studied the validity of the spin polaron picture in the frustrated J\ — J2 Heisenberg 
model. For this purpose, we have computed the hole spectral functions for the Neel, collinear, and 
disordered phases of this model, by means of the self-consistent Born approximation and Lanczos 
exact diagonalization on finite-size clusters. We have found that the spin polaron quasiparticle 
excitation is always well defined for the magnetically ordered Neel and collinear phases, even in 
the vicinity of the magnetic quantum critical points, where the local magnetization vanishes. As a 
general feature, the effect of frustration is to increase the amplitude of the multimagnon states that 
build up the spin polaron wave function, leading to the reduction of the quasiparticle coherence. 
Based on Lanczos results, we discuss the validity of the spin polaron picture in the disordered phase. 



I. INTRODUCTION 

The antiferromagnetic (AF) J\ — J2 Heisenberg Hamil- 
tonian with competing interactions to first and second 
neighbors on the square lattice represents a prototypi- 
cal model for the study of quantum phase transitions in 
low dimensional spin- 1/2 Heisenberg antiferromagnets. 
The zero point quantum fluctuations, inherent of the AF 
Heisenberg model, can be enhanced by magnetic frustra- 
tion, leading to an interesting ground state phase dia- 
gram which, after many years of investigations, has been 
quite well establishecU^i^ and it is depicted in Fig. 
[TJ As the frustration J2/J1 increases, the local magne- 
tization of the Neel state is reduced, persisting until the 
quantum critical point at about J2/J1 ~ 0.4 5 . 
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FIG. 1: Schematic magnetic phase diagram of the Ji — J2 
model. The solid lines indicates the values of the local mag- 
netizations. 

On the other hand, for J2 larger than J\, the quantum 
collinear state - neighboring spins aligned ferromagneti- 
cally along x direction and antifcrromagnctically along y 
direction, or vice versa - is stabilized by an order from 
disorder phenomenon-' 7 . On this side of the diagram, 
a decrease of J%j J\ implies an increase of frustration 
and the zero point quantum fluctuations reduce the lo- 
cal magnetization of the collinear state which vanishes at 
the quantum critical point J2/J1 ~ 0.6. In the interme- 
diate region 0.4 < Jij J\ % 0.6, quantum fluctuations are 
strong enough to destroy the magnetization, and the na- 



ture of the ground state is still under debate^, although 
the most sophisticated numerical method indicates that 
the plaquette resonance valence bound (RVB) would be 
the possible ground stated. While the experimental re- 
alization of the quantum Neel order is widely found in 
many transition metal oxides, it was only recently that 
the quasi-2D vanadium oxide Li^VOSiO^ was found as 
the first experimental realization of the J± — J2 model, 
being located at the collinear side of the phase diagram 9 . 

A little explored issue in the literature^ is the dy- 
namics of holes injected in the magnetic ground states 
of the Ji — J2 model. Since the hole strongly couples 
with the magnetic excitations, dictated by the different 
symmetries of the magnetic ground states, a strong de- 
pendence of the one-particle properties with frustration is 
expected. Hence, it is interesting to investigate whether 
the conventional quasiparticle excitations^, the spin po- 
laron picture, is valid or not in magnetically frustrated 
states. 

In the present article, we present a systematic study 
of the dependence of the one-hole spectral functions with 
frustration for the whole range of parameters. Using the 
self consistent Born approximation (SCBA) and exact di- 
agonalization, our main objective is to investigate the va- 
lidity of the spin polaron picture when the quantum criti- 
cal points are approached from the Neel and the collinear 
sides of the phase diagram. In previous investigations we 
have already studied the hole motion in the triangular 12 
and canted^ antiferromagnets. We have shown that the 
non-collincarity of the magnetic ground states, a 120° 
and tilted 180° Neel order, respectively, leads to two dif- 
ferent mechanisms for hole motion -one magnon-assisted 
and other free-like - which interfere in a non trivial way, 
giving rise to the vanishing of the quasiparticle weight 
in some regions of the Brillouin zone (BZ) 12 i 13 . In the 
present study we have found that the spin polaron pic- 
ture, which seems to be the correct description for the 
unfrustrated Neel case, remains valid for the frustrated 
Neel and collinear phases, even in the vicinity of the 
magnetic critical points. As a general feature, the effect 
of frustration is to increase the contribution of multi- 
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magnon states in the spin polaron wave function leading 
to the reduction of the QP weight, which always remains 
finite. For the Neel case, the quasiparticle dispersion 
has always the same structure, with a ground state mo- 
mentum at k = (7r/2,7r/2), even when the number of 
magnons increases as the disordered phase is reached. 
On the other hand, for the collinear phase, the structure 
of QP dispersion changes notably with frustration due to 
the interplay between the magnon assisted and free-like 
mechanisms for hole motion, leading to the divergence of 
the QP effective mass for some momenta. In the disor- 
dered regime, we have performed an exact diagonaliza- 
tion study on finite-size clusters. Guided by the SCBA, 
the exact spectra can be interpreted as the result of the 
coupling of the hole with a superposition of Neel and 
collinear short range correlations. However, it is difficult 
to extrapolate this result in the disordered phase since it 
is believed that, in the thermodynamic limit, the nature 
of the low lying magnetic excitations above the candidate 
ground stated and their eventual coupling with the hole 
could change dramatically. 

The paper is organized as follows. In Sec. II we present 
the self-consistent Born approximation, and we compare 
its predictions with exact calculations on finite-size clus- 
ters. In Sec. Ill we analyze the hole spectral functions 
obtained for the different ground states of the J\ — J2 
antifcrromagnct. Finally, in Sec. IV we state the conclu- 
sions and discuss the validity of the spin polaron picture 
in the general context of frustrated antiferromagnets. 

II. MODEL AND METHODS 

We consider the t— J\ — J 2 model to describe the motion 
of a hole in the frustrated J\ — J2 antiferromagnet, which 
is written as 

H = "* E (cU>+H.c) (1) 

<ij>cr 

+ J\ ^ Si.Sj + J 2 ^ Sj.Sfc, 

<ij> <ik> 

where < ij > and < ik > stand for summation over first- 
and second neighbors on a N sites square lattice, respec- 
tively. For simplicity, we consider only hopping terms t 
to first neighbors. The electronic operators are the pro- 
jected ones, Ci a = (1 — ni_ CT )cj CT , that obey the no double 
occupancy constraint. To compute the hole spectral func- 
tion of the t — J 1 — J 2 model in the thermodynamic limit, 
we appeal to the spinless fermion representation and the 
self-consistent Born approximatio n 14 ! 15 . 

A. Self-consistent Born approximation 

Within the SCBA scheme, the hole is represented by 
a spinless fermion hi, which carries the charge degree of 
freedom, while the magnetic excitations are treated in 



the spin- wave approximation (SWA), by means of the 
Holstein-Primakov bosons a%. Since the SWA is valid 
only for magnetically ordered states, with magnons as 
their elementary excitations, we restrict the use of the 
SCBA to the Neel and collinear phases of the J\ — J2 
antiferromagnet. Following a standard procedure^, we 
obtain the effective Hamiltonian for the hole motion 

H eff = E £ k^k^k + W 1 a q a 1 ( 2 ) 
k q 

+ ]T(A/kXV<_ k + H.c.) , 

k,q 

where hk is the Fourier transform of the spinless fermion 
operator, and a q are the Bogoliubov operators that di- 
agonalize the Heisenberg part in the SWA. 

The effective Hamiltonian comprises three terms: i) a 
free-like hopping term that takes into account the pos- 
sibility of the hole to move without disturbing the un- 
derlying magnetic background. It is characterized by the 
hole dispersion 

ffk = > ts COS — — COS K.O, 

s 2 

where ^2 S is the summation over neighbors d of a given 
site, connected by the hopping term ts (first neighbors in 
our case). Q is the magnetic wave vector, (tt,w) for the 
Neel phase, and (n, 0) or (0, ir) for the collinear phase. 
It should be noticed that a dispersive bare hole band 
will exist only for magnetic backgrounds with some non- 
collinear or ferromagnetic correlations between nearest 
neighbors, like in the collinear phases. Instead, in the 
Neel phase, e-^ — 0, and this fact will affect considerably 
the structure of the hole spectral function. 

ii) a free magnon term, characterized by the dispersion 
in the SWA 

where 

7k = \ Y Js cos2 ~T~ C0S k '^' r = \ E Js cos Q-* 5 ' 

1 S 1 1 8 

and 

lv- r . 2 Qi . . 

1 6 1 

iii) a hole-magnon interaction term, which incorporates 
the magnon-assisted mechanism for the hole motion, 
characterized by the hole-magnon vertex interaction 

M kq = -= (/3 q Vk-q - /?kUq-k) , 
V JV 

where 

/?k = tg sin — — sin k.<5, 
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and the Bogoliubov coefficients are 



' 7q - r + UJq 



2uj 




To study the hole dynamics, we compute the hole spec- 
tral function, that is, A^(uj) = — (1 / Tr)ImG^.(u>) , where 



G£(w) =< AF\h k [l/(u + 



irp 



H eff )]hl\AF > is the 



retarded hole Green function, and \AF > is the undoped 
magnetic ground state in the SWA. In the SCBA, the 
hole self-energy is given by the following self-consistent 
equation 



E k (w) = 



q 



- irf 



|Mk q | 2 ^ 

— CJ q — £k-q — S k _ q (a> — CJfc-q) 



This equation has been numerically solved for square lat- 
tices with up to 3600 sites, with a linear mesh of 10000 
frequences, and r\ + = O.Olt. Once the hole self-energy 
is obtained, the QP spectral weight can be calculated 
as Zk = [1 — dT,k(ui) /duj]~ 1 \E k , where the QP energy is 
given by the equation = £k + Ek(£'k)- 



B. Exact and SCBA comparison 
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FIG. 2: Hole spectral functions Ak(w) for the Neel state, 
N = 16, Ji/t = 0.4, k = (7r/2,7r/2), and two different values 
of the frustration, J2/J1 = 0.0 (left) and J2/J1 = 0.4 (right). 
The solid and dashed lines correspond to the Lanczos and 
SCBA results, respectively. The Lorentzian broadening for 
Lanczos curves is rj + = O.lt 

We benchmark the accuracy of the analytical SCBA by 
making a direct comparison of its predictions with Lanc- 
zos calculations performed in finite-size clusters. In Fig. 
[5] we show the Lanczos and SCBA hole spectral func- 
tions for a 4 x 4 lattice at the ground state momentum 
k = (it/2, 7t/2). The left panel corresponds to the unfrus- 
trated case, while the right panel to the frustrated case, 
J2/J1 = 0.4. It is well known in the literatur o 14 i 15 the re- 
markable agreement between SCBA and Lanczos results 



TABLE I: Exact structure factor Sq for the J\ — J2 Heisenberg 
model, with J2/J1 = 0.7, N = 16, and different values of 
the spatial anisotropy in the exchange interactions to first 
neighbors, Ji x = J\, 



k 


J\y J\X 


Jly = 0.99J1, 


Jly — 0.2Jl x 


(7T,0) 


0.81590 


1.24823 


1.57717 


(0,7T) 


0.81590 


0.38094 


0.18420 


(7T, 7r) 


0.12557 


0.13307 


0.15190 


(7r/2,7T/2) 


0.16624 


0.16635 


0.16667 


(7T,7T/2) 


0.27813 


0.35594 


0.28422 


(7T/2,7r) 


0.27813 


0.19768 


0.16669 



for the unfrustrated case, in the strong coupling regime. 
This has been ascribed to the vanishing of the first-order 
vertex corrections to the SCBA hole self- energy^. From 
Fig. [21 it is evident that the good agreement remains 
valid in the frustrated case where, despite the increase 
of zero point quantum fluctuations, the SCBA repro- 
duces the main features of the spectral function. Wc 
have found that SCBA yields quantitative good results 
for all momenta in the BZ and in the strong coupling 
limit (Ji/i < 1). 
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FIG. 3: Hole spectral function Ak(ui) for the collinear state, 
N = 16, ,h/t = 0.4, J2/J1 = 0.7, and momenta (tt,0) (up- 
per panels) and (7r/2, 7r/2) (lower panels). Left: J\ y — Ji x , 
middle: Ji y = 0.99Ji^, right: 3\ y — O.SJi^; ,h x = J\. The 
solid and dashed lines correspond to the Lanczos and SCBA 
results, respectively. 

The collinear phases break the Heisenberg Hamiltonian 
SU(2) spin-rotation symmetry, like in the Neel phase, 
and, in addition, the invariance under tt/2 real space 
rotation^. This discrete two-fold degeneracy cannot be 
spontaneously broken in Lanczos exact calculations on 
finite-size clusters, resulting in the equal presence of 
(tt, 0) and (0, 7r) correlations in the magnetic background. 
On the other hand, in SCBA this discrete degeneracy is 
broken from the onset, with the selection of one of the two 
possible magnetic wave vectors, (7r,0) or (0, 71"). There- 
fore, one should not naively compare Lanczos and SCBA 
results, until the degeneracy is somehow broken in Lane- 
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zos calculations, selecting one of the collinear orders. To 
this end, we consider spatial anisotropy in the exchange 
interactions to first neighbors, that is, different values of 
Ji along the x direction (J\ x ) and the y direction (Ji y ). 
Clearly, if J\ x is larger than J\ y , (7r,0) correlations will 
be favored. In fact, from Table I it results that even a 
slight anisotropy, J% y = 0.99Ji2;, breaks the two-fold de- 
generacy in the collinear regime and considerably favors 
(-7T, 0) correlations. 

Fig. displays the Lanczos and SCBA spectral func- 
tions, for N = 16 sites, in the collinear regime J2/J1 — 
0.7. In the isotropic case and for k = (71", 0) (upper 
left panel), there is a strong disagreement between both 
curves; the exact spectral function exhibits a two-peak 
structure but the SCBA function has only a pronounced 
peak at high energy. This disagreement can be un- 
derstood if we consider that, in each method, the hole 
moves on different magnetic background. While the 
Q = (it, 0) collinear phase has been chosen from the onset 
in the SCBA, in the Lanczos calculations (tt, 0) and (0, tt) 
correlations are present on equal footing (see Table I). 
With the introduction of a slight amount of anisotropy, 
Ji y = 0.99Jia;, the agreement between both methods 
improves considerably (middle upper panel). The small 
anisotropy induces almost no changes in the SCBA be- 
cause its magnetic background remains the same. In 
Lanczos calculations, instead, the anisotropy produces a 
noticeable qualitative change, since now the (tt, 0) corre- 
lations are largely favored over the (0, tt) ones (see Table 
I), resembling the SCBA results. For larger anisotropy, 
J\ y = 0.8 J\ x (upper right panel), the agreement between 
exact and SCBA spectral functions improves consider- 
ably, because the exact magnetic background is closer 
to the Q = (tt, 0) state assumed in the analytical calcu- 
lations. On the other hand, for momenta along the BZ 
diagonal, like (tt/2, ir/2) shown in the lower panels of Fig. 
[31 the exact hole spectral functions do not depend on the 
anisotropy, always showing a good agreement with the 
SCBA spectra. 

The performed comparison lends reliability to the 
SCBA predictions for both frustrated phases, Neel and 
collinear. In the next section, we present SCBA results 
for the hole dynamics in the thermodynamic limit. 



III. RESULTS 

A. Frustrated Neel case 

In Fig. Q] we show the SCBA spectra for the un- 
frustrated case, J2/J1 — 0, along the momentum line 
(7r/2,7r/2) -> (7r,7r) and J x /t = 0.2. For this strong 
coupling regime, we observe a strong k— dependence of 
the spectra due to the underlying magnetic Neel state. 
As a general feature, at low energy there is a spin- 
polaron QP excitation, and several resonances above it, 
the so called string excitations, whose energies scale as 
E string ~ (Ji/t) 2 ^ 3 , for all momenta. This structure of 



the spectra has been widely confirmed by several numer- 
ical and analytical technique a 15 ' 16 i 17 . While the QP ex- 
citation results from the coherent coupling of the hole 
motion with the AF spin fluctuations, the higher en- 
ergy string excitations correspond to the hole motion 
inside a linear potential generated by the overturned 
spins. The existence of a coherent QP is assigned to the 
zero point quantum fluctuations that repair pairs of mis- 
aligned spins at a characteristic time 1/ J\ , while the hole 
overturns spins at a characteristic time l/t. Although 
this polaron picture is valid for the weak coupling regime 
J > i, the very existence of the spin polaron excitations 
in the strong coupling regime t > J gives support to the 
idea that the coherence between hole motion and spin 
fluctuations is partially preserved. 




k=(Tt,Tt) 



4 k=(7t/2,Tc/2) 



FIG. 4: Hole spectral function A^(ui) for the unfrustrated 
J2/J1 ~ Neel state and Ji/t — 0.2 along the momentum 
line (7r/2,7r/2) -» (7r,7r). 

The effect of frustration is shown in Fig. [5]for J2/J1 = 
0.3. We observe a spectral weight reduction of both, the 
spin polaron and the string excitations, although the QP 
weight of the spin polaron remains always finite for the 
whole BZ, even for the extreme frustrated case (see be- 
low). Furthermore, there is a redistribution of the spec- 
tral weight at low energies which is related to the fact 
that the only mechanism for hole motion is magnon as- 
sisted. At the magnetic level, frustration induces a re- 
duction of the magnon dispersion bandwidth which can 
be associated with a reduced effective exchange J e ff- 
Therefore, the loss of QP coherence can be ascribed to 
an increase of the characteristic time of the spin fluctua- 
tions 1/J e //. In addition, it is also known that the spin- 
stiffness is smaller for the frustrated case, making the 
Neel order less rigid. Such a softening of the magnetic 
correlations explains the weakening of the linear poten- 
tial seen by the hole with the subsequent broadening of 
the string excitations. A similar effect of frustration on 
the QP excitation has been observed by Shibata et alr^ 



5 




k=(7C/2,7t/2) 



FIG. 5: Hole spectral function A^(uj) for the frustrated Neel 
state J2/J1 = 0.3 and Ji/t = 0.2 along the momentum line 

(7r/2,7r/2) -» (7T,7T). 



although their calculation, motivated by the description 
of the photoemission spectra of Sr-zCuOzCh, was per- 
formed on an extended t — t' — t" — J model. 



k=(Jt/2,7t/2) 

k=(7t,0) 

k=(jr,Ji) 




FIG. 6: Quasiparticle weight of the spin polaron as a function 
of frustration J2/J1. On the right of the vertical line a short 
range Neel state has been assumed in the calculation. 



lations with gapped magnetic excitations. The resulting 
QP weight is shown on the right of the vertical line of 
Fig. [6] In particular, we confirmed that the SCBA spec- 
tra for long and short range correlations are practically 
the same what means that hole dynamics only depends 
on the short range magnetic correlations. To understand 
the behavior of the QP weight with frustration it is in- 
structive to write down the spin polaron wave function 

„„19 



|* k >=a^ 0) ^|AF> 



q' 



( 3 ) 

where the first term represents a state with one hole and 
no magnons, the second term one hole with one magnon, 
and so on. The coefficients are the probability ampli- 
tude to have n magnons in the spin polaron state. As the 
effect of frustration is to reduce the magnon dispersion 
bandwidth, it is easier for the hole to emit and absorb 
magnons. It is then obvious that, from renormalization 
condition < Vl/kl^k >= 1, the increasing contribution of 
multimagnon states will be accompanied by the reduc- 



tion of the QP weight a£ 



1/2 



Alternatively, we can 



extend the approximated expression for the QP weight, 
derived by Kane et alM- , to the frustrated case, 



Zlc < 1 



t 2 m*z k 



J ly /1-2J 2 /Ji 



where m* is the spin polaron effective mass, to make 
evident the reduction of the QP weight with increasing 
frustration. 
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The other issue we are concerned with, is the validity 
of the spin polaron picture near the magnetic quantum 
critical point J2/J1 ~ 0.4^. In Fig. [6] we show the be- 
havior of the QP weight Zk as a function of frustration 
Ji/Jx- It is observed a monotonic decrease of the QP 
weight which remains finite for all frustration and for the 
whole BZ. Since near the critical point the local magne- 
tization is m ~ 0, the use of the spin wave theory for the 
calculation of the QP weight can be questioned. In or- 
der to have reliable results around the critical point, we 
have implemented the SCBA with a modified spin wave 
theory which only incorporates short range Neel corre- 



FIG. 7: Quasiparticle energy dispersion along the momentum 
line (0,0) — (7r,7r) — (t, 0) — (0,0), for several values of frus- 
tration in the Neel state. J2/J1 — (solid line), J2/J1 = 0.2 
(dotted line), J2/J1 = 0.4 (dashed line). Inset: bandwidth 
dispersion of the quasiparticle W (solid line) and magnetic 
W m (dashed line) excitations as a function of frustration. 

Fig. [7] displays the spin polaron dispersion for several 
values of frustration. It can be seen that the reduction 
of the QP bandwidth W is related to the reduction of 
the magnetic bandwidth W m = 2(J 2 — Ji), see inset of 
Fig. [71 This is due to the fact that the only mechanism 
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for hole motion, in the Neel state, is the magnon assisted 
one. Another consequence is the unchanged structure 
of the QP dispersion with frustration, which is also in 
accordance with the spectral weight redistribution at low 
energies stated above. The unchanged structure of the 
QP dispersion along with the reduction of the bandwidth 
and the QP weight, are all consistent with the intuitive 
picture of considering the effect of frustration on the hole 
dynamics as a loss of QP coherence due to a reduced 
effective exchange J e //- 

B. Frustrated Collinear case 

As pointed out in Section II, the hole dynamics in the 
collinear phase is quite different to the Neel phase. The 
arrangement of ferromagnetic chains perpendicular to the 
AF chains leads to an effective Hamiltonian that involves 
the coexistence of the magnon assisted and free-like pro- 
cesses for hole motion (see Eq. [5]) . The presence of such 




FIG. 8: Hole spectral function Ak(oj) for the collinear state 
Q = (n, 0), Jij J\ — 1, and Ji/t = 0.2 along the momentum 
line (ir/2,w) -> (tt/2,0). 

processes can be identified in the hole spectral functions. 
In Fig. [8] we show the SCBA prediction for the collinear 
phase Q = (7r,0) at J2/J1 = 1 and J\/t = 0.4, along the 
momentum line (ir/2, ir) — ► (tt/2, 0). For this strong cou- 
pling regime, there is again a strong k— dependence of the 
spectra due to the underlying magnetic collinear state. 
The low energy part is related to the spin polaron pic- 
ture where the hole is mainly coupled to the AF magnons, 
while at high energy there is a resonance in the spectra 
that is related to the hole motion along the ferromagnetic 
chains whose dispersion is of the tight binding form. Al- 
though the QP excitation is originated by the AF fluctu- 
ations, the string excitations practically disappear owing 
to the fact that the only retraceable paths generating a 
linear potential are one dimensional in the collinear state 




k=(7t/2,0) 



FIG. 9: Hole spectral function A-^(uj) for the more frustrated 
J2/J1 = 0.6 collinear state Q = (tt, 0) and Ji/t — 0.2 along 
the momentum line (7r/2, 7r) — > (7r/2,0). 
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FIG. 10: Quasiparticle weight of the spin polaron as a func- 
tion of frustration J2/J1. 

tically do not change when a short range collinear order 
with gapped excitations is assumed in the calculation 
near the critical point. Similarly to the Neel case, the 



. On the other hand, the 2D AF character of the Neel 
state increases noticeably the number of such retraceable 
paths, making the string excitations observable. Besides 
the reduction of the magnetic bandwidth, another effect 
of frustration in the collinear state is to tune the inter- 
k= (7l/2,0) ference between both processes for hole motion, resulting 
in a spectral weight transfer from low to high energies. 
See for instance, in Fig. [9l the merging shoulders at 
lj ~ 3t of the last three spectra, near k = (tt/2, 0), for the 
more frustrated case J2/J1 — 0.6. Although not easily 
discernible, the spin polaron excitation is well defined. 
In particular, Fig. [10] shows the monotonic reduction 
of the QP weight as the critical point J2/J1 ~ 0.51 is 
approachec— -. We again confirm that the spectra prac- 
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spin polaron picture is valid for the collinear regime, be- 
ing the QP weight always finite for the whole BZ. We 
again interpret the reduction of Zk as a consequence of 
the increasing contribution of multimagnon terms in the 
spin polaron wave function of Eq. [3] For increasing frus- 
tration -smaller J2/J1- there is a linear reduction of the 
QP dispersion bandwidth, which follows the magnetic 
one W m = 2 J 2 + J\ (see inset of Fig. [TT]), and a notable 
change of the structure of the QP dispersion (see Fig. 
[TT]). In particular, around J2/J1 ~ 2 the QP minimum 
moves from (ir,ir) and (0, tt), to (w/2,tt) leading to the 
divergence of the effective mass for some momenta. This 



1.5 



1 



1 

2 


1 1 1 








W/t 




i 


/ ' \\ 

/ ' A 

A' 1 

- / 1 ■' ' \ 




// \ \ 


5 1 1,5 


A x \ 

A x \ 
/' ■ • x \ 
//•'_*. * \ _ 


/ ■' ' *i 

■ ' V* ■ 
1 \ 


1 JI^r^L. 1/ 


/>*/' \ *. \ 



(0,0) (7t/2,7t/2) (71,71) (lt/2,lt) (0,7t) (0,0) 



k 

FIG. 11: Quasiparticle energy dispersion along the momen- 
tum line (0,0) — (tt, tt) — (tt, 0) — (0,0), for several values of 
frustration in the collinear state. J2/J1 = 2.5 (solid line), 
J2/J1 = 2 (dashed line), J2/J1 = 0.9 (dotted line), and 
J2/J1 = 0.7 (dash-dotted). Inset: bandwidth dispersion of 
the quasiparticle W (solid line) and magnetic W m (dashed 
line) excitations as a function of frustration. 

behavior is in contrast to that found in the Neel state, 
where the only effect of frustration was a redistribution 
of spectral weight at low energies. In the collinear state, 
instead, the spectral weight transfer is between the low 
and the high energy parts of the spectra which, at the mi- 
croscopic level, can be ascribed to the subtle interference 
between the magnon assisted and free like processes for 
hole motion. It is by this mechanism, ultimately tuned by 
frustration, that the QP dispersion changes its structure. 
It is interesting to note that the same behavior - namely, 
the change on the structure of the QP dispersion along 
with the non-vanishing of QP weight in the vicinity of 
the quantum critical point - has been recently observed 
in the bilayer Heisenberg antiferromagnet, where the in- 
tcrlayer coupling J± induces a quantum disordered phase 
transition^. 



C. Disordered phase 

As pointed out in the introduction, the nature of the 
disordered magnetic ground state in the intermediate 
regime 0.4 < J%j J\ ^ 0.6 is still under debate, although 
there is numerical evidence of the presence of plaque- 



tte RVB correlations^. While in the ordered regimes 
the spin wave excitations are the proper magnetic de- 
grees of freedom that couple with the hole, the lack of a 
correct analytical description of the magnetic excitations 
above the disordered phase does not allow one to envis- 
age the proper effective Hamiltonian for the hole propa- 
gation. Then, the hole dynamics in the disordered phase 
is not analytically tractable within our scheme of calcula- 
tion. Consequently, we resort to the Lanczos method to 
study the hole spectral functions on finite-size clusters. 
For these small systems, there is an important reduction 
of the disordered window relative to the thermodynamic 
limit. In particular, for N = 20 the weakening of the Neel 
and the collinear correlations occurs in a narrow window 
around J2/J1 ~ 0.58, with an enhancement of transla- 
tional symmetry breaking correlations, such as columnar 
and plaquette-like^i. In Fig. [12] we show the hole spec- 




FIG. 12: Exact hole spectral function A^(lu) for a cluster of 
TV = 20 sites evaluated at k = (tt, tt) for the Neel (dotted line), 
collinear (dashed line), and disordered (solid line) regimes. 

tral functions for N = 20 sites, momentum k = (77, it), 
and three values of frustration ranging from the Neel 
(J2/J1 = 0.57) to the collinear (J2/J1 = 0.6), across the 
disordered regime (J2/J1 = 0.58). We observe a notice- 
able spectral weight transfer from high to low energies de- 
spite the very small change in frustration, AJ2 = 0.03Ji. 
This behavior can be interpreted by recalling the differ- 
ent mechanisms for hole motion in the Neel and collinear 
regimes. The structure of the hole spectral function for 
J2/J1 — 0.57 (0.6) closely resembles that of the magnet- 
ically ordered Neel (collinear) phase, in the thermody- 
namic limit. Surprisingly, in the disordered regime, the 
hole spectral function (solid line in Fig. [T2"]) seems to be 
the superposition of the Neel and collinear cases. In par- 
ticular, as frustration increases the position of the peaks 
remains unaltered while a spectral weight transfer takes 
place. Extra features in the hole spectral functions, that 
could indicate the coupling of the hole to translational 
symmetry breaking correlations, are not observed in the 
disordered phase. It is worth noticing, however, that 
our results correspond to a small cluster size, N = 20, 
and therefore they cannot be naively extrapolated to the 
thermodynamic limit. Hence, the hole dynamics in the 
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disordered phase of the J1—J2 Heisenberg model, and the 
validity of the spin polaron picture in particular, remains 
an open problem. 

IV. CONCLUDING REMARKS 

In this work we have studied the dynamics of a hole 
injected in the magnetic ground states of the frustrated 
J\ — J2 Heisenberg model. To accomplish this task we 
have used the self consistent Born approximation and ex- 
act diagonalization on finite clusters. Our main result is 
that the spin polaron scenario remains valid whenever the 
ground state is magnetically ordered, even in the vicinity 
of the Neel and collinear quantum critical points. That 
is, there is always a well defined quasiparticle excitation 
associated with a hole dressed by magnons. Furthermore, 
the narrowing of the magnon dispersion, induced by frus- 
tration, leads to a loss of QP coherence due to the prolif- 
eration of multimagnon processes that build up the spin 
polaron wave function. 

At this point, it is interesting to discuss the physical 
scenario behind the spin polaron picture in the general 
context of 2D antiferromagnets. The main idea is based 
on the dominant pole approximation^ where the exis- 
tence of a QP pole crucially depends on the density of 
states of the low lying energy excitations which can cou- 
ple to the hole. For the simpler Neel case, there are few 
low lying spin excitations, and there is very little phase 
space available for the hole to scatter, so the low en- 
ergy states of the hole can have a long lifetime. In this 
case the underlying microscopic process for hole motion 
is the magnon assisted one. For the collinear state, the 
presence of ferromagnetic correlations induces free-like 
hole hopping processes that interfere with the magnon 
assisted ones. However the QP excitations are still well 
defined even near the magnetic critical point. For non 
collinear cases, like the triangular and canted antiferro- 
magnets, we have recently shown that the effect of such 
interference can be dramatic. In particular, for the tri- 
angular AF with positive t we have founcU^ a vanishing 
of the QP excitation in an ample region of the Brilouin 
zone. For canted Neel phases, we were able to inves- 
tigate such interference by continually tuning the cant- 
ing angle, finding a critical value at which the QP ex- 
citations disappear—. Therefore, from the present and 
previous studies, we argue that the collinearity of the 
Neel and the collinear phases are crucial for the exis- 



tence of QP excitations. This idea is also supported by 
results on the honeycomb^ 2 , and bilayer square lattices 2 ^. 
One interesting question is the breakdown of the spin 
polaron excitation and the possible realization of spin- 
charge separation^. It was suggested that this scenario 
could take place in disordered magnetic states, although 
the lack of reliable methods to compute hole spectral 
functions in the thermodynamic limit makes this task 
quite difficult. Here the nature of the low lying magnetic 
excitations is crucial. For gapped magnetic excitations, 
like in ladder— and checkerboard 2 ^ systems, a coherent 
QP motion of an injected hole is expected, at least in 
part of the BZ. On the other hand, for the kagome lattice 
there is a quite small triplet gap A ~ 0.05 J and an expo- 
nentially large number of singlet excitations within such a 
triplet gap^l. These unconventional low lying excitations 
may prevent the appearance of spin polaron excitations, 
as recently observed in numerical investigations 2 ^. In 
the present work we have numerically obtained the hole 
spectral functions for the disordered phase of the J\ — J2 
model, for N = 20 sites. We have found that these spec- 
tral functions resemble a superposition of the Neel and 
collinear cases, suggesting the existence of quasiparticle. 
Although our results cannot be naively extrapolated to 
the thermodynamic limit, we believe that if the ground 
state of the disordered regime is a plaquette RVB with 
triplet gapped excitations, the spin polaron picture could 
still be valid. 

From the experimental point of view, the quasi-2D 
vanadium oxide Li^V OSiO^ is considered a good real- 
ization of the J\ — J2 model, located at the collinear side 
of the diagram-. Recent NMR experiments 2 ^ revealed 
a very low spin dynamics which has been interpreted as 
the coexistence of (n, 0) and (0, tt) magnetic domains of 
collinear phases in the regime TV < T < J2 + J\ . Beyond 
the theoretical interest in the validity of the spin polaron 
picture, the interplay we have established between the 
microscopic mechanisms and the features of the spectral 
functions can help to interpret experimental studies 29 . 
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